Randmoness and Step-like Distribution of Pile Heights in Avalanche Models 
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The paper develops one-parametric family of the sand-piles dealing with the grains' local losses 
on the fixed amount. The family exhibits the crossover between the models with deterministic and 
stochastic relaxation. The mean height of the pile is destined to describe the crossover. The height's 
densities corresponding to the models with relaxation of the both types tend one to another as the 
parameter increases. The densities follow a step-like behaviour in contrast to the peaked shape found 
in the models with the local loss of the grains down to the fixed level [S. Liibeck, Phys. Rev. E, 62, 
(2000), 6149]. A spectral approach based on the long-run properties of the pile height considers the 
models with deterministic and random relaxation more accurately and distinguishes the both cases 
up to admissible parameter values. 

PACS numbers: 05.70Jk, 05.45Pq, 05.65+b 



I. INTRODUCTION 

In 1987 Bak ct al (BTW) introduced their sand-pile 
model [j. The model determines a system containing 
some physical quantity called sand in the original paper. 
The system is slowly loaded. Extra loading results in 
a local relaxation. The local relaxation releases energy 
that can instantly spread out to the large distances. The 
spreading mechanism is fully deterministic. The model's 
system achieves its critical state without tuning any pa- 
rameter Q. 

Numerous power laws describe the critical state. They 
have been established theoretically |3J, |4J and numerically 
0- IH 01 The model laws find their application for such 
different fields as neural networks |y, earthquakes 0, 
and solar flares [Io| . 

A great demand for the model has resulted in its mod- 
ifications. The closest versions to the original sand-pile 
are, probably, Manna's and Zhang's models 0, Hi^ . 
Manna has defined the spreading of the local relaxation 
in a stochastic way. Zhang has introduced a continuous 
sand-pile. The critical behaviour of these models exhibits 
a certain similarity. Since minor changes in the model 
rules weakly influence the critical behaviour El] an d the 
number of the changes is inexhaustible, the models need 
a strict classification. 

Many papers assign arbitrary two models to the same 
universality class if they have the same set of the ex- 
ponents determining the critical power laws 13]. The 
"heat" discussion El El El El EH based on the differ- 
ent approaches has turned to the conclusion that BTW's 
and Manna's sand-piles belong to the different universal- 
ity classes EH- Preliminary investigation suggests that 



'Finance Academy under the Government of the Russian Federa- 
tion.; Electronic address: shapoval@mccme.ru 

tlnstitut de Physique du Globe de Paris; Electronic address: 
shnir@mitp.ru, shnirman@ipgp.jussieu.fr 



BTW's and Zhan g's s and-piles should represent the same 
universality class [lflj . 

The paper [2(| has introduced the family of the mod- 
els realizing the crossover between Manna's and Zhang's 
sand-piles. The control parameter deals with the energy 
of the local relaxation. The small values of the energy 
characterized Manna model while the extremely big val- 
ues lead to Zhang- type behaviour. The properties of the 
sand distribution over the system reflect the crossover. 

The local relaxation is defined for j2fl] 'sand-piles as the 
loss of the energy down to the fixed level. The family of 
the sand-piles in [2l| deals with the loss of energy on the 
fixed amount. 

According to (2l| , its family of the models makes the 
crossover between BTW's sand-pile and the random walk 
crossing some modification of Manna's sand-pile. This 
crossover corresponds to the relatively small energy men- 
tioned above. On the other hand, the family determines 
a sophisticated limit behaviour as the energy tends to 
infinity. The classification based on the power laws fails 
to describe a great diversity appeared in this continuous 
sand-pile family. The paper 22] has introduced an ap- 
propriate global functional, whose evolution calculated in 
terms of the spectrum determines the system dynamics. 

The paper [2j| introduces the energy propagation 
mechanism that lacks the local symmetry, thus leading to 
the sand-piles with the quenched disorder. These models 
depend on some parameter exhibiting the value of the 
asymmetry. BTW sand-pile belongs to the family of the 
models and corresponds to the absense of the asymme- 
try. Establishing the deterministic relaxation the l2^l's 
family of the models principally differs from that of |2df 
and [21j . The properties of this family are found but the 
models' crossover to BTW sand-pile is not considered. 

In this paper we investigate the sand-pile family of |2l| . 
The family's sand distribution is found to be principally 
different from that appearing in the usually discussed 
models, in particular, in Zhang's and Ref. [20j]'s models. 
As the control parameter is relatively big the difference 
between the deterministic and stochastic relaxation al- 
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most disappears. Then the sand distribution over the 
system is proved to become similar for the both types of 
the models. However the spectral properties select the 
sand-piles with the deterministic relaxation. 



II. MODEL 



The model deals with a two-dimensional square lattice 
L x L. Each cell contains hij grains, where hu is less 
than an integer threshold H . At each time moment a 
cell is chosen at random. Its number of the grains 
(further, height) hij increases on 1: 



hi 



1. 



If the resulting height hij remains less than H, then noth- 
ing more happens at the moment. Otherwise, the cell 
becomes unstable and relaxes. Relaxation depends 
on an integer control parameter n. The unstable cell dis- 
tributes its n grains "in equal parts" among its 4 nearest 
neighbours. Namely, as n — 4fc each neighbour gets ex- 
actly k grains. 

Naturally, there exist the numbers n = 4k + r, where 
the residue r < 4 is not equal to zero. Then each neigh- 
bour gets k grains and the rest r grains are distributed 
to four different neighbours at random. The following 
formula expresses the idea of the construction: 



^ncighbour(i, j) 
^noighbour(i, j) 



"'neighbour^, j) ' 
^ncighbour(i, j) ~t~ k 



1. 



During this relaxation other cells can achieve the 
threshold H and become unstable. Then they relax ac- 
cording to the same rules. If a boundary cell relaxes, 
[n/4] or [n/4] + 1 grains leave the lattice and dissipate, 
where [x] is the integer part of x. (The dissipation is 
bigger for the corner cells). 

The sequential acts of relaxation are called an avala- 
nche. The size of any avalanche is the number of the 
unstable cells during the avalanche counting according 
to their multiplicity. The dissipation at the boundary 
assures that the avalanches are defined correctly and their 
size is finite. 

The case of n = H = 4 corresponds to the original 
sand-pile of . It worth noting that the paper [2(| intro- 
duces relaxation with the loss of any unstable cell down 



to the fixed level. Namely, hi 



H 



with this hi 



grains passed one-by-one to the nearest neighbours; the 
receiver for each grain is found at random. These changes 
in the rules qualitatively influence the systemThe limit 
behaviour (n — > oo) of 20] 's family and the models in 
question is different. 

The family really depends only on parameter n. The 
values of H does not influence the system dynamics. 
They determine the admissible interval [H — n, H) of the 
heights. 
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FIG. 1: Normalized density of the heights; the dashed lines 
are the linear fits, L — 256. 



III. PILE HEIGHT DENSITIES 

Following the ideas of [2(j we establish some features of 
the sand distribution over the lattice. Let the normalized 
heights hij n be hij n = (hij—H+n)/n. Then hij n G [0, 1). 
Let a function p n (-) be the density of the normalized 
heights hij n : Pn(k/n) is the number of /iy n being equal to 
k/n. Then the normalized densities p n (-)n follow four flat 
steps with a high accuracy (fig. OJ for n = 4k. The steps 
correspond to the values of the density p4 for BTW's 
sand-pile. (The points of p4 shown in fig. ^ are in good 
agreement with the analytical values found in |24|). 

Each density is fitted by a linear function (fig.HJ. The 
slopes slightly increase as n goes up and must be satu- 
rating as n tends to infinity. 

The following construction manages to compare quan- 
titatively four steps of p„(-) as n > 4 with the four values 
of pa{-) representing BTW sand-pile. Given n = 4, 5, . . ., 
the values of the density p n (-) is sampled into four 
bins [0,0.25), [0.25,0.5), [0.5,0.75), and [0.75,1) with re- 
ported values tpn(-) at 0, 0.25, 0.5, and 0.75 respectively. 
For example, ip n {0) = p„(0) + p n {l/n) + . . . + p n (jo/n), 
where jo is the biggest integer with jo/n < 0.25. In par- 
ticular, coincides with pi{-). Then i/} n ( m ) is defined 
in 4 points. Each value ip{k/4), k = 0, 1, 2,3, represents 
one step of /?„(•). 

Let 



k=0 




Then a n measures the difference between the steps and 
the function pi(-) corresponding to BTW sandpilc. 

The difference from BTW sand-pile increases as a n 
goes up (fig. [5J • The reported values of <j n indicate that 
there exist a certain maximum of the difference corre- 
sponding to n ~ 60. This observation agrees with the 
change of the tendency exhibited by the sand-pile family 
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FIG. 2: Difference between the steps of p n and the function 
pi measured by the functional tr n . 
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FIG. 3: Normalized density of the heights; the dashed lines 
are the linear fits, L — 256. 



of 20] for the intermediate values of the control param- 
eter. 

In the same way, the normalized densities are intro- 
duced to develop the models with n =/= Ak starting with 
n = 3. The computer experiment proves (fig-EJ that the 
densities p n quickly (as n goes up) achieves four steps of 
P4. The linear fits have slopes appearing close to that for 
Pik- 

The following reasoning gives a rough explanation of 
the step-like behaviour of the densities. Given the pa- 
rameter value n, the heights are naturally divided into 
four intervals [H — n,H — 3n/4), [H - 3n/4, H — n/2), 
[H — n/2, H—n/4), [H — n/4, H) since each act of relax- 
ation maps one interval into another (possibly, except- 
ing the boundary values due to the random effect). The 
normalization (on n) transforms these intervals to the 
domain of definition of four steps in fig-^ and|3] 

It worth noting that the sand-piles of [12, l20j demon- 
strate peaked densities of the average height in contrast 
to our p n . 



So, in terms of the sand distribution the model family 
exhibits a certain similarity. There exist possibilities to 
select the limit behaviour (n — > 00) and define an inter- 
mediate deviation from the extreme cases. The densities 
p n for n — Ak and n ^ Ak become almost indistinguish- 
able as n has the order of dozens. 



IV. SPECTRAL APPROACH 

Another approach gives evidence of the diversity of the 
two cases (n = Ak and n 7^ Ak). It deals with the spec- 
trum of the average height h = L~ 2 . —1 hij. The av- 
erage height is calculated at the end of each time moment 
and treated as the function on time, h(t). The paper 
has successfully used h(t)'s spectrum to describe the pa- 
per's sand-pile's dynamics. Our /i(t)'s spectrum appears 
to be noisy, therefore it is averaged over its several real- 
izations. Besides, each h(t)'s realization is stored in the 
bins of some length A. Then the Fourier transform de- 
termines the spectrum £. The highest frequency is 1/A 
in the case. 

A formal procedure applied for the spectrum calcula- 
tion consists of four steps. 

1. For some fixed N the time moments, which h(t) 
is catalogued at, is divided into iV non- intersecting 
intervals of the same length T. 

2. Given any fixed interval of step 1, the averaging 
of /i(f)-values over relatively small sub-intervals is 
applied to reduce the number of data for further 
numerical application of the fast Fourier transform. 
Let A be the length of the small sub-intervals and 
r be their quantity, rA = T . Then a signal Xj is 
defined as the arithmetic mean of all the /i(£)-values 
in the j-th small sub-interval, j — 1, 2, . . . , r. 

3. The Fourier transform determines the spectrum of 
x-j. Let 



v 3=1 



t))t 



where 



is the mean of ; 



Xj, J 



= 1, 



Then the frequencies fk are fc/(rA), k = 
and the power spectrum of the signal Xj is 

defined as £(/ fc ) = \x k \ 2 . 

A. Each of iV intervals defined in step 1 generates its 
own spectrum £(fk)- The averaging of £(/fc) results 
in the stabilized spectrum, which the construction 
aims at. The stabilized spectrum depends on the 
model parameter n. Thus, the notation £ n (fk) is 
keeping further for the stabilized spectrum. 

According to [2l| , the family of the spectra admits a 
certain normalization. If A is proportional to n, the exhi- 
bition of £n(fk)/ n versus dimensionless frequencies /^A 
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FIG. 4: Normalized spectrum £ n /n vs normalized frequencies 
f = fk- The inset contains the graphs in the boxed part of 
the figure; A = 8n, N = 16, r = 8192, L = 256. 



collapses the graphs on the interval of the high frequen- 
cies and remains reasonable on the other interval (fig.0J. 

Leaving the main part of the spectrum as it is we stress 
attention on some interval of the low frequencies. The 
corresponding part of the graphs is boxed in fig. 0] and 
put in its inset. In this part the graphs (for n = 4k) 
establish a noticeable rise as n goes up to the value being 
close to 60. Then the tendency principally changes. The 
graphs turn down to the spectrum £4 at their left part 
while demonstrating a definite peak (fig. 0] inset). It 
agrees with the behaviour of a n reported earlier. 

The spectrum's propagation to the left needs much 
bigger h(t)'s domain of definition that is hardly possi- 
ble. On the other hand, the highest achieved frequencies 
correspond to the avalanches of the biggest size observed 
during computer simulation. Bigger avalanches are prac- 
tically not observable even for much longer simulation. 
Then the computer experiment must be producing the 
constant spectrum just at the left of the obtained graphs' 
part in fig. 0] 

The horizontal coordinates of the investigated box are 
not absolute constants. They depend on the lattice 
length L and the simulation time. As L increases the 
box moves to the left becoming invisible in fig. 0] 



The influence of the simulation time is essential. The 
investigated part of the spectrum reflects the frequen- 
cies of the rare, big, and strongly dissipative avalanches. 
The lack of the data limits the results covering these 
avalanches. It concerns the partial conclusions about the 
rare avalanches in BTW's and Manna's sand-pile |25j . 

In contrast to the quick convergence of p n n, n ^ 4k, to 
their four values, the family £ n /n of the spectra exhibits 
quite a different tendency. Since the high-frequency spec- 
tra are rather similar (fig. 03a), the analysis is focused on 
the low frequencies (fig. EJb) corresponding to the boxed 
part in fig. [SJa. 

As n = 3 the normalized spectrum has its horizontal 
interval. While n increasing, £ n /n changes on this inter- 
val. Fig. demonstrates the minor changes as n = 30 
and completely another behaviour as n is equal to 150 
and 301. So, for the simulated big ns the low- frequency 
spectrum essentially deviates from £3 as well as from £4^ 
for big k. 

According to fig. 01 and [3] the long-time evolution of 
the sand-pile exhibits a great complexity. In contrast to 
the space features, several patterns do not exhaust the 
spectrum behaviour. Up to the developed experiments 
the sand-piles with the deterministic and stochastic re- 
laxation have the different spectra. 



V. CONCLUSION 

Summarizing, we develop the family of the sand-piles. 
The control parameter n is the number of the grains that 
any unstable cell passes to its neighbours. As n = 4k, the 
propagation of the grains through the lattice is fully de- 
terministic despite the models with n 7^ Ak involve some 
random effect. In terms of the sand distribution the ran- 
dom effect disappears as n is sufficiently big. However 
the trace of the deterministic relaxation remains visible 
in terms of the average height's spectrum, selecting some 
peak at low-frequency spectrum. With their peculiar evo- 
lutionary properties the sand-piles corresponding to the 
deterministic relaxation may admit a certain prediction. 
This hypothesis agrees with the effective precursors found 
for BTW's sand-pile in g^. 
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